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ABSTRACT: Weak gravitational lensing provides a unique method to map directly the dis- 
tribution of dark matter in the universe and to measure cosmological parameters. This cosmic- 
shear technique is based on the measurement of the weak distortions that lensing induces in the 
shape of background galaxies as photons travel through large-scale structures. This technique 
is now widely used to measure the mass distribution of galaxy clusters and has recently been 
' detected in random regions of the sky. In this review, we present the theory and observational 

status of cosmic shear. We describe the principles of weak lensing and the predictions for the 
shear statistics in favored cosmological models. Next, we review the current measurements 
of cosmic shear and show how they constrain cosmological parameters. We then describe the 
prospects offered by upcoming and future cosmic-shear surveys as well as the technical challenges 
Cf*^ . that have to be met for the promises of cosmic shear to be fully realized. 

6 ■ 1 INTRODUCTION 

Gravitational lensing provides a unique method to directly map the distribution 
of dark matter in the universe. It relies on the measurement of the distortions 
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that lensing induces in the images of background galaxies (for reviews, see Bartel- 
^ ■ mann &; Schneider 1999, Bernardeau 1999, Kaiser 1999, Mellier 1999, Narayan & 

Bartelmann 1996, Schneider 1995, Wittman 2002). This method is now widely 
used to map the mass of clusters of galaxies (see Fort & Mellier 1994 for a review) 
and has been extended to the study of superclusters (Gray et al. 2002, Kaiser et 
al. 1998) and groups (Hoekstra et al. 2001). Recently, weak lensing was statisti- 
cally detected for the first time in random patches of the sky (Brown et al. 2003; 
Bacon et al. 2002; Bacon, Refregier & Ellis 2000; Hamana et al. 2003; Hammerle 
et al. 2002; Hoekstra et al. 2002a; Hoekstra, Yee & Gladders 2002a; Jarvis et 
al. 2002; Kaiser, Wilson & Luppino 2000; Maoli et al. 2001; Refregier, Rhodes, 
& Groth 2002; Rhodes, Refregier & Groth 2001; van Waerbeke et al. 2000; van 
Waerbeke et al. 2001a; van Waerbeke et al. 2002; Wittman et al. 2000). These 
cosmic shear surveys provide direct measurements of large-scale structure in the 
universe and, therefore, of the distribution of dark matter. Unlike other meth- 
ods that probe the distribution of light, weak lensing measures the mass and 
can thus be directly compared to reliable theoretical models of structure forma- 
tion. Cosmic shear can therefore be used to measure cosmological parameters 
in the context of these models, thereby opening wide prospects for cosmology 
(Bernardeau, van Waerbeke & Mellier 1997; Hu & Tegmark 1999; Jain & Seljak 
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Figure 1: Illustration of the effect of weak lensing by large-scale structure. The 
photon trajectories from distant galaxies (right) to the observer (left) are deflected 
by intervening large-scale structure (center). This results in coherent distortions 
in the observed shapes of the galaxies. These distortions, or shears, are on the 
order of a few percent in amplitude and can be measured to yield a direct map 
of the distribution of mass in the universe. 

1997; Kaiser 1998). 

In the present review, we describe the theoretical and observational status of 
cosmic shear. Earlier reviews of this fast-evolving field can be found in Hoekstra 
et al. (2002b), Mellier et al. (2001), van Waerbeke et al. (2002), Wittman 
(2002). Here, we first describe the principles of weak lensing (Section 2). We then 
summarize the different statistics used to measure cosmic shear and describe how 
they are used to constrain cosmological parameters (Section 3). In Section 4, we 
survey the different methods used to derive the lensing shear from the shapes 
of background galaxies. We present, in Section 5, the current observations and 
their cosmological significance. Future cosmic-shear surveys and the prospects 
they offer for cosmology are described in Section 6. In Section 7, we outline how 
systematic effects present challenges that must be met for the potential of these 
future surveys to be fully realized. We summarize our conclusions in Section 8. 

2 THEORY 

The idea of cosmic shear can be traced back to a lecture given by Richard Feyn- 
man at Caltech in 1964 (J.E. Gunn, personal communication). Several theo- 
rists (e.g., Gunn 1967, Jaroszyhsky et al. 1990, Kristian & Sachs 1966, Lee h 
Paczyhsky 1990, Schneider & Weiss 1988) then studied the propagation of light 
in an inhomogeneous universe. Predictions for the statistics of the weak-lensing 
distortions were then computed in a modern cosmological context by several 
groups (Babul & Lee 1991, Blandford et al. 1991, Kaiser 1992, Miralda-Escude 
1991, Villumsen 1996). More recently, the power of cosmic shear to measure 
cosmological parameters was the object of many theoretical studies (Bernardeau, 
van Waerbeke & Mellier 1997; Jain & Seljak 1997; Hu & Tegmark 1999; Kaiser 
1998; Kamionkowski et al. 1997; van Waerbeke, Bernardeau & Mellier 1999). In 
this section, we briefly describe the principles of weak lensing and show how this 
technique can be used to map the dark matter in the universe. 

As they travel from a background galaxy to the observer, photons get deflected 
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Figure 2: Illustration of the geometrical meaning of the shear 7^ and of the 
ellipticity e^. A positive (negative) shear component 71 corresponds to an elon- 
gation (compression) along the x-axis. A positive (negative) value of the shear 
component 72 corresponds to an elongation (compression) along the x = y axis. 
The ellipticity of an object is defined to vanish if the object is circular (center). 
The ellipticity components ei and €2 correspond to compression and elongations 
similar to those for the shear components. 

by mass fluctuations along the line of sight (see Figure 1). As a result, the appar- 
ent images of background galaxies are subject to a distortion that is characterized 
by the distortion matrix: 



where 59i(9) is the deflection vector produced by lensing on the sky. The conver- 
gence k is proportional to the projected mass along the line of sight and describes 
overall dilations and contractions. The shear 71 (72) describes stretches and com- 
pressions along (at 45° from) the x-axis. Figure 2 illustrates the geometrical 
meaning of the two shear components. 

The distortion matrix is directly related to the matter density fluctuations 
along the line of sight by 



where $ is the Newtonian potential, x is the comoving distance, Xh is the co- 
moving distance to the horizon, and <9j is the comoving derivative perpendicular 
to the line of sight. The radial weight function g(x) is given by 



where r = ci~ 1 Da, and Da is the angular-diameter distance. The function n(x) 
is the distribution of the galaxies as a function of the comoving distance x from 
the observer and is assumed to be normalized as / dxn(x) = 1. 

As we discuss in Section 4, galaxy shapes can be averaged over a patch of the sky 
to measure the shear, which is thus an observable. The shear pattern expected in 
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Figure 3: Shear map derived by ray-tracing simulations by Jain, Seljak & White 
(2000). The size and direction of each line gives the amplitude and position angle 
of the shear at this location on the sky. The displayed region is 1° x 1° for an 
SCDM (Einstein-De Sitter) model. Tangential patterns about the overdensities 
corresponding to clusters and groups of galaxies are apparent. A more complex 
network of patterns is also visible outside of these structures. The root- mean- 
square shear is approximately 2% in this map. (From Jain et al. 2000) 

a standard Cold Dark Matter (SCDM) model is shown in Figure 3 for a 1 x 1 deg 2 
region. Jain, Seljak, & White (2000) derived this map from ray tracing through 
N-body simulations. Tangential patterns around the overdensities corresponding 
to clusters and groups of galaxies, along with a more complicated network of shear 
fluctuations, are apparent. By inverting the lensing equation (Equation 2), the 
shear map can be converted into a map of the projected mass k and, therefore, 
of the dark matter distribution. 

3 COSMIC-SHEAR STATISTICS AND COSMOLOGY 

The statistical characteristics of the cosmic-shear field can be quantified using a 
variety of measures, which can then be used to constrain cosmological models. 
First, we consider the most basic two-point statistic of the shear field, namely the 
two-dimensional power spectrum (Jain & Seljak 1997, Kaiser 1998, Kamionkowski 
et al. 1997, Schneider et al. 1998a). The shear power spectrum C\ is defined as a 
function of multipole moment I (or inverse angular scale) by Sf=i (^(1)^(1')) = 
(27r) 2 (5(l — l')C\, where tildes denote Fourier transforms (with the conventions of 
Bacon, Refregier, & Ellis 2000), S is the two-dimensional Dirac-delta function, and 
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Figure 4: Shear power spectrum for different cosmological models and for source 
galaxies at z s = 1. The SCDM model is COBE normalized and thus has a higher 
amplitude than the three cluster-normalized models ACDM, OCDM, and rCDM. 
The thin dashed line shows the ACDM spectrum for linear evolution of structures. 
Notice that for I > 1000 (corresponding approximately to angular scales 8 < 10') 
the lensing power spectrum is dominated by nonlinear structures. 

the brackets denote an ensemble average. Applying Limber's equation in Fourier 
space (e.g., Kaiser 1998) to Equation 2 and using the Poisson equation, one can 
easily express the shear power spectrum C\ in terms of the three-dimensional 
power spectrum P(k,x) of the mass fluctuations dp/p and obtain 




where a is the expansion parameter, and H$ and £l m are the present value of 
the Hubble constant and matter density parameter, respectively. The lensing 
power spectra for four CDM models are shown in Figure 4 (see color insert) 
(see Bacon, Refregier &: Ellis 2000 for the exact cosmological parameter values 
of each model). They were derived using the fitting formula for the nonlinear 
matter power spectrum P(k,x) of Peacock & Dodds (1996). In Figure 4, the 
galaxies were assumed to lie on a sheet at a redshift z s = 1. A more realistic 
redshift distribution n(z) would require corrections of only approximately 10% 
on these power spectra, as long as the median redshift of the galaxies were kept 
at z s = 1. 

The three cluster-normalized models (ACDM, OCDM and rCDM) yield power 
spectra of similar amplitudes but with different shapes. The COBE-normalized 
SCDM model has more power on a small scale and thus yields a larger normal- 
ization. For the ACDM model, the power spectrum corresponding to a linear 
evolution of structures is also shown in Figure 4 for comparison. For / > 1000 
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(corresponding to angular scales smaller than approximately 10'), nonlinear cor- 
rections dominate the power spectrum (Jain & Seljak 1997), making cosmic-shear 
sensitive to gravitational instability processes. 

The measurement of the lensing power spectrum can thus be used to measure 
cosmological parameters, such as Q m , J7a, 08, and T. (Bernardeau, van Waerbeke 
& Mellier 1997; Hu & Tegmark 1999; Jain & Seljak 1997; Kaiser 1998; van Waer- 
beke, Bernardeau & Mellier 1999). A full-sky cosmic-shear survey would yield 
a precision of these parameters comparable to that for future cosmic microwave 
background (CMB) missions. More realistically for the short term, a precision 
on the order of 10% can be achieved with surveys of approximately 10 square 
degrees. Such cosmic-shear surveys can also be combined with CMB anisotropy 
measurements to break degeneracies present when the CMB is considered alone 
(Hu & Tegmark 1999; Contaldi, Hoekstra k Lewis 2003). This would yield im- 
provements in the precision of cosmological parameters by approximately one 
order of magnitude. In addition, the use of photometric redshifts can provide 
a tomographic measurement of matter fluctuations and improve the precision of 
cosmological parameters by up to an order of magnitude (Hu 1999, Hu & Keeton 
2002, Taylor 2001). 

In practice, it often is more convenient to measure other two-point statistics. 
In particular, the variance of the shear a 2 = (j 2 ) in randomly placed cells is 
widely used. It is related to the shear power spectrum by 



where W\ is the Fourier transform of the cell aperture. The shear two-point cor- 
relation functions (eg. Kaiser 1998) are Fourier transforms of the power spectra 
and have also been measured by various groups (eg. Bacon et al. 2002; van 
Waerbeke et al. 2001a; Hoekstra et al. 2002a). The M ap statistics (Schneider et 
al. 1998a) is another convenient statistic based on the average of the shear within 
a compensated filter. Its advantages are that its window function in /-space is 
narrow, that adjacent cells are effectively uncorrelated, and that it can be easily 
related to the statistics of the projected mass k. 

In analogy with electromagnetism, a tensor field like the shear field 7j can 
be decomposed into an electric (E, or gradient) component and a magnetic (B, 
or curl) component (Crittenden et al. 2001b, Kaiser 1992, Kamionkowski et al. 
1997, Stebbins 1996). Because the distortions it induces arise from a scalar field 
(the gravitational potential), weak lensing only produces E-type fluctuations. On 
the other hand, systematic effects and intrinsic galaxy alignments are likely to 
produce both E-type and B-type fluctuations (See Heavens 2001 for a review; 
also see Section 7.6). The presence of B-modes can thus be used as a measure of 
contaminants to the cosmic-shear signal. In practice, this decomposition can be 
performed using the M ap statistic (Schneider et al. 1998a, van Waerbeke et al. 
2001a) or the power spectrum (Hu & White 2000, Padmanabhan, Seljak & Pen 
2002, Pen et al. 2001). 

As is apparent in Figure 3, the shear field is not Gaussian but, instead, has 
apparent coherent structures. Bernardeau, van Waerbeke, & Mellier (1997) have 
shown that the skewness 

«? - (K3) (6) 
S *-J^2 (6) 

of the convergence field k can be used to break the degeneracy between ag and 
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Cl m , which is present when the shear variance alone is considered. The measure- 
ment of higher-point statistics is thus of great cosmological interest, but their 
computation is made difficult by the fact that the angular scales accessible to 
obervation {0 < 10') are in the nonlinear regime. Predictions for higher-point cor- 
relation functions have been calculated using perturbation theory, the hierarchi- 
cal ansatz, (Bernardeau & Valageas 2000; Bernardeau, van Waerbeke & Mellier 
1997; Hui 1999; Munshi & Coles 2000; Munshi & Jain 2001; van Waerbeke et 
al. 2001b) and halo models (Cooray & Hu 2001a; Cooray, Hu, & Miralda-Escude 
2000). These techniques have also been used to compute the full probability dis- 
tribution function of the convergence field (Bernardeau & Valageas 2000, Munshi 
Sz Jain 2000, Valageas 2000) and the errors in two-point statistics (Cooray &: Hu 
2001b; Munshi & Coles 2002; Schneider et al. 2002). 

Another possible way to compute the full nonlinear field is to use ray tracing 
through N-body simulations (Blandford et al. 1991; Hamana, Colombi & Mellier 
2000; Jain, Seljak & White 2000; Premadi et al. 2001; Wambsganss, Cen, & 
Ostriker 1998; White & Hu 2000) to produce simulated shear maps such as those 
by Jain, Seljak & White (2000) (shown in Figure 3). These can be used to 
compute and study other proposed measures of non-Gaussianity such as peak 
statistics (Jain & van Waerbeke 2000), a generalized maximum likelihood (Taylor 
& Watts 2000), and cluster counts (Bartelmann, King & Schneider 2001). In 
practice, the complex geometry of surveys makes it difficult to infer a convergence 
k map from the observed shear 73. For this reason, a number of researchers have 
recently proposed the use of high-order statistics of the shear field rather than 
the skewness S3 of the convergence (Bernardeau, van Waerbeke, & Mellier 2003; 
Schneider & Lombardi 2002; Takada & Jain 2002; Zaldarriaga & Scoccimarro 
2002). 

4 SHEAR MEASUREMENT METHODS 

Because the sought-after lensing signal is of only a few percent in amplitude, the 
data acquisition and analysis must be performed carefully, and systematic effects 
must be tightly controlled. All the current measurements of cosmic shear were 
derived from deep optical images taken with charged-coupled devices (CCD). It 
is advantageous for the exposures to be homogeneous in depth and for the ground 
to be subject to as small a seeing as possible. The fields are generally chosen to 
lie far away from each other to ensure that they are statistically independent and 
to minimize cosmic variance. 

The first step in the data analysis is image processing. After flat fielding, 
the different exposures are co-added to produce the final reduced images. If 
necessary, any instrumental distortion induced by the telescope optics is corrected 
for at this stage. This can be done very accurately by measuring the astrometric 
offsets from several dithered exposures. For a detailed description of the different 
image-processing steps, see Kaiser et al. (1999). An example of a processed deep 
image from the cosmic-shear survey of Bacon, Refregier & Ellis (2000) is shown 
in Figure 5. 

The next step consists of deriving an estimator for the shear from the shapes 
of the galaxies in the co-added images. The point-spread function (PSF), which 
smears the images of galaxies and is generally not circular, complicates this task. 
In general, the PSF varies spatially and in time, and it must be measured and 
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Figure 5: Example of an deep image in the cosmic-shear survey by Bacon, Re- 
fregier & Ellis (2000). This corresponds to a 1 h exposure with the EEV camera 
on the William Herschel Telescope (WHT). The field of view is 8' x 16' and 
achieves a magnitude depth of R ~ 26 (5a detection). The bright objects are 
saturated stars. The faint objects comprise approximately 200 stars and approxi- 
mately 2000 galaxies that are usable for the weak-lensing analysis. (From Bacon, 
Refregier & Ellis 2000) 
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modelled for each image individually. This can be done by measuring the shape 
of the stars in the field, whose number can be optimized by tuning the galactic 
latitudes of the observations. Figure 6 shows the ellipticity pattern of the PSF 
from one of the William Hershel telescope (WHT) fields of Bacon, Refregier & 
Ellis 2000. 

Several methods have been developed to tackle this difficult and crucial task. 
The more rigorous method of Kaiser, Squires & Broadhurst (KSB; 1995), further 
developed by Luppino & Kaiser (1997) and Hoekstra et al. (1998), replaced the 
earlier method by Bonnet & Mellier (1995). The KSB method is now widely used 
for cluster studies, and it has been used by the majority of the groups involved in 
measuring cosmic shear (see Bartelmann &: Schneider 1999 for a detailed review 
of the KSB method). It is based on the measurement of the quadrupole moment 
of the galaxy surface brightness /(x), 

Qij = J (fxxiXjw(x)I(x.), (7) 

where w(x) is a weight function conveniently taken to be Gaussian. These mo- 
ments capture the lowest-order shape information of the galaxy and can be com- 
bined to form the ellipticity of the galaxy 

e _ Qu — Q22 e _ 2Q12 ^ 
Qu + Q22 Qu + Q22 

The first (second) component of the ellipticity describes compressions and elon- 
gations along (at 45° from) the x and y axes (see Figure 2 for an illustration). 
The ellipticity vanishes for circular galaxies. 

The first step in the KSB method is to correct the observed galaxy ellipticity 
ef for the anisotropy of the PSF. The corrected galaxy ellipticity e 9 is given by 

^ ^ ~ ^sm(-fsm) e ) (9) 

where e* is the PSF ellipticity derived from the stars, and Pg m and P* m are the 
smear susceptibility tensors for the galaxy and star, respectively, and can be 
derived from higher moments of the images. The shear in a patch of the sky can 
then be measured by averaging over the (corrected) ellipticities in the patch of 
the sky using 

7 = (p 7 rv>, (10) 

where the tensor P 7 quantifies the susceptibility to shear acting before isotropic 
PSF smearing and is given by 

p^pi-pKPLr'pL, (ii) 

where the shear susceptibility tensors P^ h and P s * h for the galaxies and the stars 
can also be measured from higher moments of their respective light distribution. 

The KSB method was thoroughly tested using realistic simulated images by 
Erben et al. (2001) and Bacon et al. (2001a). These studies showed that it is 
accurate to within a few tenths of percent in reconstructing an input shear and 
is thus sufficient for the current cosmic-shear surveys (see Section 5). However, 
this accuracy is insufficient for future, more sensitive surveys (see Section 6). 
Moreover, Kuijken (1999) and Kaiser (2000) have shown that the KSB method 
is ill-defined mathematically and unstable for PSF's found in practice. 
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Figure 6: PSF ellipticity pattern measured in WHT field shown in the Figure 5: 
(a) PSF ellipticities measured from the stars in Figure 5. (b) Ellipticity residuals 
for those stars after the KSB corrections. In this survey, the root-mean-square 
stellar ellipticity is found to be approximately 7% before the corrections and 
negligible after the correction. (From Bacon, Refregier Sz Ellis 2000) 



This inadequacy has led a number of researchers to develop alternative meth- 
ods. Rhodes, Refregier & Groth (2000) have modified the KSB method to be 
better suited for HST images. Kuijken (1999) considered a different approach 
that consisted of fitting the observed galaxy shape with a smeared and sheared 
circular model. Kaiser (2000) introduced a new method based on a finite resolu- 
tion shear operator. Refregier & Bacon (2003) and, independently, Bernstein & 
Jarvis (2001) developed a new method based on the decomposition of the galaxies 
into shape components or "shapelets" . The gauss-hermite orthogonal basis func- 
tions used in this approach allow shears and PSF convolutions to be described as 
simple matrix operations, using the formalism developped for the quantum har- 
monic oscillator (Refregier 2003). These new methods are promising but require 
extensive testing to establish whether they will achieve the required precision. 

5 OBSERVATIONS 

Because the expected distortions are only of a few percent, the measurement 
of cosmic shear requires large survey areas and excellent image quality. Early 
searches for cosmic shear signals with photographic plates were unsuccesful (Kris- 
tian 1967; Valdes, Jarvis & Tyson 1983). Mould et al. (1994), performed the 
first attempt to detect a cosmic-shear signal with CCDs, but only derived an 
upper limit. Using the same data, Villumsen (1995) reported a 4.5<r detection. 
Schneider et al. (1998b) then reported a detection of cosmic shear in one of three 
QSO fields, an area too small to draw any constraints on cosmology. 

Within a few weeks, four independent groups (Bacon, Refregier & Ellis 2000; 
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Kaiser, Wilson & Luppino 2000; van Waerbeke et al. 2000; Wittman et al. 2000) 
reported the first firm statistical detections of cosmic shear using three different 
4m-class telescopes: the Cerro Tololo Inter- American Observatory (CTIO), the 
Canada- France-Hawaii Telescope (CFHT), and the William Herschel Telescope 
(WHT). These were later confirmed by more precise measurements of the cosmic- 
shear amplitude from the ground (Bacon et al. 2002; Brown et al. 2003; Maoli et 
al. 2001; Hamana et al. 2003; Hoekstra et al. 2002a; Hoekstra, Yee &: Gladders 
2002a; Jarvis et al. 2002; van Waerbeke et al. 2001; van Waerbeke et al. 2002a) 
using 2-m - (MPG/ESO), 4-m- (WHT, CFHT, CTIO) and 8-m- (Very Large 
Telescope, Keck, Subaru) class ground-based telescopes. Meanwhile, cosmic shear 
was also detected (Hammerle et al. 2002; Rhodes, Refregier & Groth 2001) and 
measured (Refregier, Rhodes & Groth 2002) from space using HST. Space-based 
surveys are currently limited by the small field of view of HST, but they are deeper 
and less prone to systematics thanks to the absence of atmospheric seeing. 

Table 1 summarizes the existing cosmic-shear surveys and highlights the wide 
range of telescopes, survey areas, and depths. The shear variance (Equation 5) 
measured recently by several groups is shown in Figure 7 (see color insert) as a 
function of the radius a circular cell. The results by Hammerle et al. (2002) 
and Hamana et al. (2002a) are not displayed. Note that, for the shear variance, 
the data points at different angular scales are not independent. For comparison, 
the shear variance predicted for a ACDM model with <rg = 1 is also shown in 
Figure 7. The model is displayed for median galaxy redshifts of z m from 0.8 to 
1.0, corresponding approximately to the range of depths of the top five surveys 
displayed (van Waerbeke et al. 2002a; Brown et al. 2003; Bacon et al. 2002, 
WHT and Keck; Refregier et al. 2002). These observations are approximately 
consistent with each other and with the ACDM model on angular scales from 0.7 
to 20 arcmins. This is compelling given that these were performed with different 
telescopes (and therefore different instrumental systematics) and independent 
data-analysis pipelines. 

The bottom two surveys (Hoekstra et al. 2002b; Jarvis et al. 2002) have a 
median redshift in the range z m ~ 0.6-0.7 and yield lower shear variances. As 
indicated by the theoretical curves in figure 7, this redshift dependence is expected 
in CDM models and thus confirms the detection of a cosmic shear signal. 

Although the shear variance was displayed for the purpose of Figure 7, some 
of the groups have used instead the shear correlation function or M ap statistic 
to quantify their lensing signal. Recently, several groups used the latter statistic 
to separate their signal into E and B components (see Section 3) and thus to 
estimate and subtract systematic effects (Jarvis et al. 2002; Hamana et al. 2003; 
Hoekstra, Yee & Gladders 2002; van Waerbeke et al. 2001, 2002a). Pen, van 
Waerbeke & Mellier (2001) and Brown et al. (2003) chose instead to measure 
directly the shear power spectrum for each E and B mode. Figure 8a shows 
the measurement of the E and B signals using the M ap statistic by Hoekstra, 
Yee &: Gladders (2002a), for several i? c -magnitude ranges in the Red-Sequence 
Cluster (RCS) survey. They measure a clear lensing signal apparent as E modes. 
However, the significant B modes reveals the presence of residual systematics 
on small scales (#<10'), especially for bright galaxies. In their analysis, these 
authors use the amplitude of the B modes as a measure systematic uncertainties. 

Existing cosmic-shear measurements already yield interesting constraints on 
the amplitude of the matter power spectrum o"§ on which the lensing signal 
strongly depends. For instance, Figure 8b shows cosmological constraints for 
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Figure 7: Shear variance as a function of the radius 6 of a circular cell. The 
data points correspond to recent results from different groups: van Waerbeke et 
al. (2002a), Brown et al. (2003), Bacon et al. (2002, WHT and Keck), Refregier 
et al. (2002), Hoekstra et al. (2002b), Jarvis et al. (2002). When relevant, the 
inner error bars correspond to noise only, whereas the outer error bars correspond 
to the total error (noise + cosmic variance). The measurements by Hammerle 
et al. (2002) and Hamana et al. (2003) are not displayed. The solid curves 
show the predictions for a ACDM model with $l m = 0.3, as = 1, and T = 0.21. 
The galaxy median redshift was taken to be z m = 1.0,0.9, and 0.8, from top to 
bottom, respectively, corresponding approximately to the range of depth of the 
top five surveys. The bottom two surveys (Hoekstra et al. 2002b; Jarvis et al. 
2002) have a median redshift in the range z m ~ 0.6-0.7 and, as expected, yield 
lower shear variances. 
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Figure 8: Cosmic shear measurement by Hoekstra, Yee & Gladders (2002a) from 
the RCS survey, (a): measurement of the E (top panels) and B (bottom panels) 
modes of the M ap statistics variance as a function of aperture scale 6. Several 
ranges of the i? c -magnitude are considered. The errors correspond to la sta- 
tistical uncertainties. The presence of significant B Modes reveal the presence 
of residual systematics on small scales (9 < 10'). (b): Cosmological constraints 
derived from their measurement of the E-mode M ap statistics for 22 < R c < 24 
after accounting for the residual 5-modes. The joint constraints are shown for 
f2 m and as in the ACDM model. Priors from CMB and galaxy surveys were used 
to marginalize over f2tot> I\ and the source redshift z s . The contours indicate 
68.3%, 95.4%, and 99.9% confidence levels. (From Hoekstra, Yee & Gladders 
2002a) 
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Table 1: current cosmic-shear surveys 



Reference 


Telescope 


A / 1 V \ 

Area (deg J 


Mag. limit 


°8 a 


Wittman et al. 2000 


CTIO 


1.0 


i?<26 




van Waerbeke et al. 2000 


CFHT 


1.7 






Kaiser et al. 2000 


CI HI 


c\r' 

0.96 


T C\ A T t /or 

7<24, K <25 




oacon et ai. zuuu 


Will 


U.o 


xi ^ ZD 


L5U -0.50 


Maon et al. 2001 


VL1 


0.65 


1 < 24.5 


i no+0.036 


Rhodes et al. 2001 


TTnm /(I 71 1 I ~) /"A ,A 

Hb 1 /Wr PC2 


a Ar 

0.05 


i < 26 


n m +0.25 
°- 91 -0.30 


van Waerbeke et al. 2001a 


Ur ri 1 


0.0 


J ^ Z4.0 


n oo+0.02c 
U -°°-0.02 


Hammerle et al. 2002 


rib 1 /b lib 


n no 
0.02 




llUCAotldi ct djl. iUUia 


CFHT CTIO 


94 


R < 94 


81 + - 07 
U ' Oi -0.09 


van Waerbeke et al. 2002a 


CFHT 


8.5 


7<24.5 


o.98i8:gg 


Refregier et al. 2002 


HST/WFPC2 


0.36 


(J) ~ 23.5 


o.94i8:^ 


Bacon et al. 2002 


WHT, Keck 


1.6 


i?<26 


o.97i8:Ji 


Hoekstra et al. 2002b 


CFHT, CTIO 


53 


i?<24 


o.86l8:8 5 


Brown et al. 2003 


MPG/ESO 


1.25 


i?<25 


n 79 +0.09 
u -' z -0.09 


Hamana et al. 2003 


Subaru 


2.1 


i?<26 


69+ - 18 
U - Dy -0.13 


Jarvis et al. 2002 


CTIO 


75 


i?<23 


n 71 +0.06 
u - ' i -0.08 



"for ACDM model with Q m — 0.3, !T2a = 0.7; V is marginalised over or set to 0.21 when possible; 
errors correspond to 68%CL 

for combination of Maoli et al. (2001), van Waerbeke et al. (2000), Bacon et al. (2000), and 
Wittman et al. (2000); cosmic variance not included. 

c cosmic variance not included. 



a ACDM model derived by Hoekstra, Yee & Gladders (2002a) from the measure- 
ment of the M ap statistic for their 22 < R c < 24 galaxy sample. The degeneracy 
between <7g and O m apparent in the figure is typical of cosmic-shear measurements 
involving only two-point statistics. Hoekstra, Yee & Gladders (2002a) found that 
the constraints are well described by a 8 n ^ 2 = 0A6t° f 7 (95% CL), where priors 
from CMB and galaxy survey data have been used to marginalize over T and Qa- 
Other surveys find constraints of the same form with similar exponents for Q m . 

Table 1 lists the values of ag found by the different groups. The errors have 
been converted to 68% confidence level (CL) when necessary, and values of £l m = 
0.3, fi a = 0.7, and T = 0.21 have been assumed, when possible. The most 
recent values are displayed in Figure 9 (see color insert), along with their average 
erg ~ 0.83 ± 0.04 (lcr). The values derived by the different cosmic shear groups 
are in the range ~ 0.7-1.0, with the most recent measurements (Brown et al. 
2003; Hamana et al. 2003; Jarvis et al. 2002) yielding lower values. The 2-3cr 
mutual inconsistencies between some of the measurements may be symptomatic of 
a small level of residual systematics. In particular, calibration errors in the shear 
measurement method would not be detected via the E-B decomposition and are 
a likely explanation (see Hirata & Seljak 2003). Another source of discrepancy 
arises from the different fitting functions for the non-linear evolution of the power 
spectrum (see section 2). Apart from Brown et al. (2003; see also Contaldi et al. 
2003) who used the more accurate results of Smith et al. (2003), the other groups 
used the Peacock & Dodds (1997) fitting function which yields an underprediction 
of as by roughly 5-10%. The impact of these systematics will be discussed further 
in section 7. 
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Figure 9: Comparison of the determination of a$ by different groups and methods. 
The errors have all been converted to la, and a ACDM model with Vt m = 0.3 and 
r = 0.21 was assumed when possible. The vertical dotted lines show the average 
as — 0.83 ±0.04 of the seven cosmic-shear measurements and associated la error. 
The normalization from cluster abundance (e.g., Pierpaoli et al. 2002) as well as 
that derived from CMB anisotropies and local large-scale structure (galaxy and 
Lyman a surveys; Spergel et al. 2003) are also shown. 

Interestingly, an independent measurement of as is provided by the abun- 
dance of X-ray clusters and can be directly compared to this value. Initially, this 
method yielded normalisations of as ~ 0.9 (Eke et al. 1998; Pierpaoli, Scott & 
White 2001; Viana & Liddle 1999), in agreement with the early cosmic shear re- 
sults (see table 1 and figure 9). This value was subsequently revised downward to 
as ~ 0.7-0.8 by the use of the observed, rather than simulated, mass-temperature 
relation in clusters (Borgani et al. 2001; Reiprich & Bohringer 2001; Seljak 2001; 
Viana, Nichol & Liddle 2001; Pierpaoli, Scott k White 2002). Recently, Spergel 
et al. (2003) derived a value of as = 0.84 ± 0.04 (68%CL) from a joint analy- 
sis of CMB anisotropy measurements from the Wilkinson Microwave Anistropy 
Probe (WMAP) and other experiments, galaxy clustering and the Lyman a for- 
est. For comparison, these results along with a representative value of the revised 
cluster-abundance normalization (Pierpaoli, Scott &: White 2002) are also shown 
in figure 9. The average of the cosmic shear results is formally in good agreement 
with the determination of as using the other techniques. Future surveys are how- 
ever needed to confirm this, by resolving the discrepancies between the current 
cosmic shear measurements. Also, a full likelihood analysis would be required 
to establish the significance of the agreement between the different techniques 
(see Contaldi et al. 2003). This comparison is important as it provides a strong 
test of the ACDM model, the gravitational instability paradigm, the physics of 
clusters, and of the biased formation of galaxies. 

As discussed in Section 3, the degeneracy between as and O m can be broken 
by measuring higher-order correlation functions of the lensing field. Bernardeau, 
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Mellier, & van Waerbeke (2002) recently reported the first detection of a non- 
Gaussian shear signal using the 3-point shear correlation function formalism of 
Bernardeau, van Waerbeke & Mellier (2003). Another measure of non-gaussianity 
was recently performed by Miyazaki et al. (2002) using peak statistics in their 
Subaru survey. Although these results are consistent with that expected from 
structure formation models, larger survey areas are needed to infer cosmological 
constraints. 

Another approach to probe the dark matter is to measure the bias between 
the mass and galaxies by cross-correlating the shear map with that of the light 
from foreground galaxies in the same region of the sky (Cooray 2002, Schneider 
1998, van Waerbeke 1998). Hoekstra, Yee & Gladders (2001a), Hoekstra et al. 
(2002b), and Wilson, Kaiser & Luppino (2001) have recently measured this cross- 
correlation on large scales. 

6 FUTURE SURVEYS AND PROSPECTS 

The existing measurements described above are primarily limited by statistics. 
They will therefore be improved upon by ongoing surveys on existing telescopes, 
such as the Legacy Survey on CFHT (CFHTLS; Mellier et al. 2000), the Deep 
Lens Survey (Wittman et al. 2002), surveys with the Subaru telescope (Hiroyasu 
et al. 2001) and the Sloan Digital Sky Survey (Stebbins, McKay & Frieman 1995). 
Future instruments dedicated to surveys and for which cosmic shear is a primary 
science driver are being planned, such as Megacam on CFHT (Boulade et al. 
2000), the Visible and Infrared Survey Telescope for Astronomy (VISTA; Taylor 
et al. 2003), the Large aperture Synoptic Survey Telescope (LSST; Tyson et al. 
2002a,b), or the novel Panoramic Survey Telescope and Rapid Response System 
(Pan-STARRS; Kaiser, Tonry & Luppino 2000). From space, the new Advanced 
Camera for Surveys (ACS) on HST and, much more ambitiously, the future 
Supernova Acceleration Probe satellite (SNAP; Perlmutter et al. 2003; Rhodes 
et al. 2003; Massey et al. 2003; Refregier et al. 2003) also offer exciting prospects. 
Table 2 lists the characteristics of some of these future surveys. Broadly speaking, 
ground-based measurements will cover large areas, whereas space-based surveys 
will yield higher-resolution maps and reduced systematics thanks to the absence 
of atmospheric seeing. 

Radio surveys offer another interesting prospect. Ongoing efforts are aimed at 
detecting cosmic shear with the FIRST radio survey (Chang & Refregier 2002, 
Kamionkowski et al. 1997, Refregier et al. 1998). The future radio telescopes 
LOFAR (Low Frequency Array) and SKA (Square Kilometer Array) will yield 
cosmic-shear measurements of comparable sensitivity to the most ambitious op- 
tical surveys (Schneider 1999). The advantages of radio surveys are that they 
cover large solid angles, that the bright radio sources are at a higher redshift, 
and that the PSF is fully predictable and reproducible. 

These future surveys will provide very accurate measurements of cosmolog- 
ical parameters through the measurement of the lensing power spectrum and 
higher-order statistics. Figure 10 (see color insert) shows, for instance, the ex- 
cellent accuracy with which the lensing power spectrum will be measured with 
the SNAP wide survey. They will also allow us to test some of the foundations 
of the standard cosmological model. For example, the measurement of the power 
spectrum on nonlinear scales at different redshift slices, and of the hierarchy of 
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Table 2: Future cosmic-shear surveys 



Telescope/Survey 


Ground/Space 


Diameter 


FOV 


Area 


Start" 


Ref. b 






(m) 


(deg 2 ) 


(deg 2 ) 






DLS 


ground 


2x4 


2 x 0.3 


28 


1999 c 


1 


CFHTLS 


ground 


3.6 


1 


172 


2003 


2 


VST 


ground 


2.6 


1 


xl00 d 


2004 


3 


VISTA 6 


ground 


4 


2 


10000 


2007 


4 


Pan-STARRS 


ground 


4 x 1.8 


4x4 


31000 


2008 


5 


LSST 


ground 


8.4 


7 


30000 


2012 


6 


SNAP 


space 


2 


0.7 


300 


2011 


7 



a planned start date of the cosmic shear surveys 

6 references: 1: Wittman et al. 2002; 2: Mellier et al. 2000; 3: Kuijken et al. 2002; 4: Taylor et al. 

2003; 5: Kaiser et al. 2000; 6: Tyson et al. 2002a,b; 7: Perlmutter et al. 2003, Rhodes et al. 2003, 

Massey et al. 2003, Refregier et al. 2003 

c the survey will be complete in 2003 

d a survey of several 100 deg 2 is planned 

e assuming the availability of the optical camera 



high-order correlation functions, will yield a direct test of the gravitational in- 
stability paradigm. The lensing power spectrum can also be used to measure 
the equation of the state of the dark energy w and thus complement supernovae 
measurements in the constraining of quintessence models (Benabed & Bernardeau 
2001; Hui 1999; Huterer 2001; Hu 2001, 2002; Munshi & Wang 2002; Weinberg 
&; Kamionkowski 2002). Figure 10 shows, for instance, that a change of 40% in 
w can easily be measured by SNAP (see Refregier et al. 2003). Cosmic-shear 
measurements can also be used to test general relativity (Uzan & Bernardeau 
2001). 

Another promising approach to measure weak lensing is to use the fluctuations 
of the CMB temperature as the background sources (Bernardeau 1997, 1998; 
Cooray & Kesden 2002; Hirata & Seljak 2002; Hu 2000; Seljak 1996; Seljak & 
Zaldarriaga 1999; Zaldarriaga & Seljak 1998). Because these fluctuations are 
produced at a redshift of approximately 1100, CMB lensing provides a probe of 
the evolution of mass fluctuations at redshifts larger than those probed by optical 
galaxies. Lensing indeed produces distinct non-Gaussian signatures that can be 
used to reconstruct the foreground mass distribution and probe the growth of 
structures. This approach will become feasible with the advent of future CMB 
missions such as Planck Surveyor or ground-based instruments with high angular 
resolutions. 

7 CHALLENGE: SYSTEMATIC EFFECTS 

For the future surveys described above to yield their full potential, a number of 
challenges must be met. Indeed, these surveys require the measurement of shears 
on the order of 1% with an accuracy better than 0.1%; thus, they require a very 
tight control of systematic effects. In the following, we review the main sources 
of systematics. 
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Figure 10: Prospects for the measurement of the weak-lensing power spectrum 
with future weak-lensing surveys. The solid line represents the lensing power 
spectrum of a ACDM model with f2 m = 0.3 (and an equation of state w = — 1). 
The boxes correspond to the band-averaged la errors for the SNAP wide survey 
(300 deg 2 , 100 galaxies per arcmin with a root-mean-square intrinsic ellipticity 
dispersion of o~ e = 0.31). The precision is highest in the nonlinear region [see 
the linear ACDM power spectrum (dotted line)} and will thus provide a test of 
the gravitational instability paradigm. This model can easily be distinguished 
from models that differ from it by a 17% change in Q m (dashed line) or by a 40% 
change in the dark-energy equation of state w (dot-dashed line). 
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7. 1 PSF Anisotropy 

The most serious systematic for ground-based surveys is that produced by the 
rather large PSF ellipticities (~ 10%) observed by the different groups. The 
KSB method (see Section 4) provides a correction of the PSF anisotropy by, at 
most, a factor of approximately 10, which is insufficient for a precision of 0.1% 
in shear. Although new shear measurement methods may improve upon this (see 
Section 4 and below), the correction is fundamentally limited by the finite surface 
density and signal-to-noise ratio of the stars used to measure the PSF shape. It 
is, therefore, necessary for the PSF ellipticity to be reduced in hardware as well 
as in software. This can be done by using tighter constraints in the tracking and 
optical systems of the new telescopes and instruments as well as with interactions 
between telescope designers, engineers, and scientists. 

7.2 Shear Measurement Method 

Even if the PSF ellipticity were guaranteed to be isotropic, future measurements 
will be limited by the precision of the shear measurement methods. The KSB 
method is accurate to measure shears of approximately 1% with only a 10% 
accuracy (Bacon et al. 2001, Erben et al. 2001). This will soon become an 
important limiting factor (see Hirata & Seljak 2003). Although new methods 
(see Section 4) promise to improve upon this, more extensive simulations are 
required to establish that the same accuracy can be reached for shears of 0.1%. 

7.3 Redshift Distribution 

To convert cosmic-shear measurements into constraints on cosmological param- 
eters, the redshift distribution of the background galaxies must be known (see 
Section 4). The uncertainty in the median galaxy redshift is already one of the 
dominant contributions to the uncertainty in the amplitude of the matter power 
spectrum from cosmic shear. The determination of the galaxy redshift distribu- 
tion is made difficult by the depth of the cosmic-shear surveys. In addition, the 
sample of galaxies is not simply magnitude limited, but it is subject to complex 
selection cuts throughout the shear measurement pipeline. Dedicated spectro- 
scopic surveys and photometric redshift studies (such as that by Brown et al. 
2003) are thus required to overcome this limitation. 

7.4 CCD Nonlinearities 

All the shear measurement methods rely on the assumption that the instrumental 
response is linear. It is therefore important to test whether CCD cameras do not 
have subtle pixel-to-pixel nonlinearities that would induce biases in the shear 
measurements. The mean shear offset and the shear gradient across the chip 
found by van Waerbeke et al. (2001) could be due to this effect. 

7.5 Overlapping Isophotes 

Two neighboring galaxies in an image yield "peanut-shaped" isophotal contours 
and thus appear to have aligned ellipticities. This overlapping-ellipticities effect 
may produce a spurious ellipticity correlation signal on small scales, van Waer- 
beke et al. (2000) suggested that this effect may explain the excess shear signal, 
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which they measured on small scale and which disappeared when close pairs were 
discarded. Evidence for this effect was also found in the image simulation by Ba- 
con et al. (2001). More extensive and detailed simulations would need to be 
performed to ascertain and calibrate this effect. 

7.6 Intrinsic Correlations 

The measurement of the lensing shear relies on the assumptions that, in the ab- 
sence of lensing, the ellipticities of the galaxies are uncorrelated. However, an 
intrinsic correlation of galaxy shapes could exist owing to the coupling of the 
galaxy angular momentum or shape to the tidal field or to galaxy interactions 
(see Heavens 2001 for a review). Theoretical estimation of the size of this effect 
has been performed using numerical simulations (Croft & Metzler 2000; Heav- 
ens, Refregier &; Heymans 2000; Jing 2002) and analytical methods (Catalan, 
Kamionkowski & Blandford 2001; Crittenden et al. 2001a,b; Lee & Pen 2001; 
Mackey, White & Kamionkowski 2001). Measurements of intrinsic correlations 
have also been performed (Brown et al. 2000; Pen, Lee & Seljak 2000). Although 
considerable uncertainty remains regarding the amplitude of this effect, a consen- 
sus is arising that intrinsic correlations are likely to be small for the deep current 
surveys with z m ~ 1, but may be dominant for shallower surveys with z m <0.2 
(see however the conflicting results of Jing 2002). Although lensing distortions 
are coherent over a large redshift range, intrinsic alignments are only significant 
for small physical separations. Photometric redshifts can thus be used to separate 
and reduce intrinsic correlations from cosmic-shear signals (Heymans & Heavens 
2002; King & Schneider 2002a,b). Another approach is to search for the -B-type 
correlation signal produced by intrinsic correlations (Crittenden et al. 2001b). 

7.7 Theoretical Uncertainties 

Most of the signal in cosmic-shear surveys arises from small scales (9 < 10') and 
thus from nonlinear structures (Jain & Seljak 1997). The existing prescriptions 
for computing the nonlinear corrections (Peacock & Dodds 1997; Ma 1998) to 
the matter power spectrum are only accurate to approximately 10% and disagree 
at that level with one another (see discussion in Huterer 2001). This theoretical 
uncertainty will soon become one of the dominating errors in the determination of 
cosmological parameters from cosmic-shear surveys (see discussion in van Waer- 
beke et al. 2002). New prescriptions based on more recent N-body simulations 
such as those by Smith et al. (2003) will help improve the accuracy of the pre- 
dictions. The problem of the prediction of higher-order statistics is even more 
difficult, but it is not as pressing given the large uncertainties in the the measure- 
ments within current surveys (see Bernardeau, Mellier, van Waerbeke 2002). The 
inclusion of second-order terms in the weak-lensing approximation will eventually 
also be required (Cooray 2002, Cooray & Hu 2002). 

These systematic effects must be controlled to match the statistical accuracy 
of future surveys. Along with the studies suggested above, further measurements 
with different colors and comparison of various surveys in overlapping regions 
would help to control systematics and to test data analysis pipelines. On the 
theoretical side, larger simulations coupled with advanced analytical techniques 
will be required. 
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8 CONCLUSIONS 

Cosmic shear has emerged as a powerful method to measure the large-scale struc- 
ture in the universe. It can be thought of as the measurement of background 
fluctuations in the space-time metric. Although other methods rely on assump- 
tions relating the distribution of light to that of the mass, weak lensing is based 
on "clean" physics and can be directly compared to theory. The past three years 
have yielded impressive observational progress, as the first statistical detections 
and measurements of cosmic shear have been achieved. In analogy with the 
CMB, cosmic shear has moved from the COBE era to that of the first generation 
of anisotropy experiments. However, the measurement of cosmic shear differs 
from that of CMB anisotropics in several respects. First, the fluctuations are 
on the order of 10 -2 as opposed to 10 -5 , making them easier to measure while 
retaining the validity of linear calculations. Second, the cosmic-shear field is 
non-Gaussian and therefore contains more information than that quantified by 
the power spectrum. 

Existing cosmic-shear measurements have started to yield significant constraints 
on cosmological parameters. The measurement of the amplitude of the matter 
power spectrum o^O^ 5 from cosmic shear should soon replace that derived from 
the local abundance of clusters. This latter technique is indeed limited by the 
finite number of bright clusters and by systematic uncertainties in the physics of 
clusters. With better statistics, the angular and redshift dependence of the shear 
signal as well as with higher-order moments of the convergence field will break 
the degeneracy between eg and O m and yield constraints on further parameters 
such as Qa and V. 

The current measurements of cosmic shear are now primarily limited by statis- 
tics. They will therefore be improved upon by a number of upcoming and future 
instruments such as Megacam, VST, VISTA, LSST, and Pan-STARRS from the 
ground, and HST/ACS and SNAP in space. For several of these instruments, a 
weak-lensing survey has been listed as one of the primary science drivers. These 
surveys will potentially yield measurements of cosmological parameters that are 
comparable in precision and complementary to those derived from the CMB. 
They will also be able to address more far-reaching questions in cosmology by 
measuring parameters beyond the standard model. For instance, they can be 
used to provide a test of the gravitational instability paradigm, a measure of the 
equation of state of the dark energy, and a test of general relativity. 

For these instruments to yield the full promise of cosmic shear, a number of 
challenges have to be met. First, observationally, systematic effects, such as the 
PSF anisotropy and CCD nonlinearities, must be controlled and corrected for. 
From the theoretical point of view, calculations of the nonlinear power spectrum, 
of high-order statistics, and of the associated errors must be improved to meet the 
precision of future measurements. The observational and theoretical efforts re- 
quired to overcome these difficulties are worthwhile given the remarkable promise 
that cosmic shear offers to cosmology. 
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